The dynamical Casimir effect in braneworlds 
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In braneworld cosmology the expanding Universe is realized as a brane moving through a warped 
higher-dimensional spacetime. Like a moving mirror causes the creation of photons out of vacuum 
fluctuations, a moving brane leads to graviton production. We show that, very generically, KK- 
particles scale like stiff matter with the expansion of the Universe and can therefore not represent 
the dark matter in a warped braneworld. We present results for the production of massless and 
Kaluza-Klein (KK) gravitons for bouncing branes in five-dimensional Anti de Sitter space. We 
find that for a realistic bounce the back reaction from the generated gravitons will be most likely 
relevant. This Letter summarizes the main results and conclusions from numerical simulations which 
are presented in detail in a long paper [M. Ruser and R. Durrer, Phys. Rev. D 76, 104014 (2007), 
larXiv:0704.0790] . 
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Introduction: String theory, the most serious can- 
didate for a quantum theory of gravity, predicts the ex- 
istence of 'branes', i.e. hypersurfaces in the 10- (or 11- 
) dimensional spacetime on which ordinary matter, e.g. 
gauge particles and fermions, are confined. Gravitons 
can move freely in the 'bulk', the full higher dimensional 
spacetime [lj. 

The scenario, where our Universe moves through a five- 
dimensional Anti de Sitter (AdS) spacetime has been 
especially successful in reproducing the observed four- 
dimensional behavior of gravity. It has been shown that 
at sufficiently low energies and large scales, not only grav- 
ity on the brane looks four dimensional |2fl , but also cos- 
mological expansion can be reproduced J3| . We shall con- 
centrate here on this example and comment on behavior 
which may survive in other warped braneworlds. 

We consider the following situation: A fixed 'static 
brane' is sitting in the bulk. The 'physical brane', our 
Universe, is first moving away from the AdS Cauchy hori- 
zon, approaching the second brane. This motion corre- 
sponds to a contracting Universe. After a closest en- 
counter the physical brane turns around and moves away 
from the static brane. This motion mimics the observed 
expanding Universe. 

The moving brane acts as a time-dependent boundary for 
the 5D bulk leading to production of gravitons from vac- 
uum fluctuations in the same way a moving mirror causes 
photon creation from vacuum in dynamical cavities 
Apart from massless gravitons, braneworlds allow for a 
tower of Kaluza-Klein (KK) gravitons which appear as 
massive particles on the brane leading possibly to phe- 
nomenological consequences. 

We postulate, that high energy stringy physics will lead 
to a turnaround of the brane motion, i.e., provoke a re- 
pulsion of the physical brane from the static one. This 
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motion is modeled by a kink where the brane velocity 
changes sign. As we shall see, a perfect kink leads to 
divergent particle production due to its infinite acceler- 
ation. We therefore assume that the kink is rounded 
off at the string scale L s . Then particles with energies 
E > E s = X/L s are not generated. This setup represents 
a regular 'bouncing Universe' as, for example the 'ekpy- 
rotic Universe' Four-dimensional bouncing Universes 
have also been studied in Ref. Q. 

Moving brane in AdSs: Our starting point is the met- 
ric of AdSs in Poincare coordinates: 



ds" = gABdx dx 
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The physical brane (our Universe) is located at some 
time-dependent position y = yt(t), while the static brane 
is at a fixed position y = y s > yb(t). The scale factor on 
the brane is 
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where we have introduced the brane velocity v and the 
conformal time 77 on the brane. If v <C 1, the junction 
conditions lead to the Friedmann equations on the brane. 
For reviews see 0, Q. Defining the string and Planck 
scales by K5 = and K4 
(RS) fine tuning condition 0] implies 
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We assume that the brane energy density is dominated 
by a radiation component. The contracting (t < 0) and 
expanding (t > 0) phases are then described by 
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where H = (da / drf) / a 2 is the Hubble parameter and we 
have used that rj ^ t if v -C 1. A small velocity also 
requires yb(t) <C L. The transition from contraction to 
expansion is approximated by a kink at t — 0, such that 
at the moment of the bounce 

\v(0)\ = = £ , a b = a(0) = 4= , ff b 2 = ^| . (5) 



Tensor perturbations: We now consider tensor per- 
turbations hij on this background, 



ds 2 = — [-dt 2 + (Sij + 2h ij )dx t dx j + dy 2 ] 
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For each polarization, their amplitude h satisfies the 
Klein-Gordon equation in AdSs [8] 
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h{t,y;k)=0 
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where k = |k| is the momentum parallel to the brane and 
h is subject to the boundary (2nd junction) conditions 



(vd t +d y )h\ Vb{t) = -► dyh\ yb{ t) 



and dyh\y a 



= . 

Being interested in late-time (low energy) effects, we have 
approximated the first of those conditions by a Neumann 
condition (v <C 1)- Then, the spatial part of Eq. (J7J) 
together with (JSJ) forms a Strum-Liouville problem at any 
given time and therefore has a complete orthonormal set 
of eigenfunctions {(fi a (t, y)}^Lcr These 'instantaneous' 
mode functions are given by 



Mt) = 

<f>n{t,y) = 
C v (m,x,y) = 



y s yb(t) 



(9) 
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N n (t)y 2 C 2 (m n (t),y b (t),y) with 
Yi(mx)J v (my) — Ji(mx)Y u (my) (10) 



and satisfy [— <9 2 + (3/y)d y ](f> a (y) — m 2 t (f> a (y) as well as 
([8]). iV„ is a time-dependent normalization condition. 
More details can be found in Q. The massless mode 4>q 
represents the ordinary four-dimensional graviton on the 
brane, while the massive modes are KK gravitons. Their 
masses are quantized by the boundary condition at the 
static brane which requires Ci(m n ,yb,y s ) — 0. At late 
times and for large n the KK masses are roughly given 
by m n ~ nir jy s . The gravity wave amplitude h may now 
be decomposed as |j[ 



h{t,y;k) = J y§ ^ <? a ,k y) (H) 



where the prefactor assures that the variables q a ,k are 
canonically normalized. Their time evolution is deter- 
mined by the brane motion [cf. Eq. [|14p]. 
Localization of gravity: From the above expressions 



and using L/yb(t) = a(t), we can determine the late- 
time behavior of the mode functions <j) a on the brane 
(y b « I « !/ s ) 
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At this point we can already make two crucial ob- 
servations: First, the mass m n is a comoving mass. 
The instantaneous energy of a KK graviton is Lu n ^ — 
\/k 2 + m^, where k denotes comoving wave number. 
The 'physical mass' of a KK mode measured by an ob- 
server on the brane with cosmic time dr = adt is therefore 
m n /a, i.e. the KK masses are redshifted with the expan- 
sion of the Universe. This comes from the fact that m n 
is the wave number corresponding to the y direction with 
respect to the bulk time t which corresponds to conformal 
time r] on the brane and not to physical time. It implies 
that the energy of KK particles on a moving AdS brane 
is redshifted like that of massless particles. From this 
alone we would expect the energy density of KK modes 
on the brane decays like 1/a 4 . 

But this is not all. In contrast to the zero mode 
which behaves as 4>o{t,yb) oc 1/a the KK-modc functions 
^n(^2/fe) decay as 1/a 2 with the expansion of the Uni- 
verse and scale like 1 / ^/yl- Consequently the amplitude 
of the KK modes on the brane dilutes rapidly with the 
expansion of the Universe and is in general smaller the 
larger y s . This can be understood by studying the prob- 
ability of finding a KK-graviton at position y in the bulk 
which turns out to be much larger in regions of less warp- 
ing than in the vicinity of the physical braneQ. If KK 
gravitons are present on the brane, they escape rapidly 
into the bulk, i.e., the moving brane looses them, since 
their wave function is repulsed away from the brane. This 
causes the additional 1/a-dependence of 4> n (t,yb) com- 
pared to 4>o(t,yb). The l/.y/2/7-dependence expresses the 
fact that the larger the bulk the smaller the probabil- 
ity to find a KK-graviton at the position of the moving 
brane. This behavior reflects the localization of gravity: 
traces of the five-dimensional nature of gravity like KK 
gravitons become less and less 'visible' on the brane as 
time evolves. As a consequence, the energy density of 
KK gravitons at late times on the brane behaves as 



Pkk oc 1/a 6 



(13) 



It means that KK gravitons redshift like stiff matter and 
cannot be the dark matter in an AdS braneworld since 
their energy density does not have the required 1 /a 3 be- 
havior. They also do not behave like dark radiation 
as one might naively expect. This new result is derived 
in detail in Ref. 0. It is based on the calculation of 
(A 2 (i,?/ fc ,k)) ex (q 2 a ^{t))4> 2 a (t,y b ) where the bracket in- 
corporates a quantum expectation value with respect to 
a well-defined initial vacuum state and averaging over 
several oscillations of the field 0. An overdot denotes 
the derivative with respect to t. The scaling behaviour 
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(fT3| is due to (j>%{t,yb) only, since graviton production 
from vacuum fluctuations has ceased at late times (like 
in radiation domination) which is necessary for a mean- 
ingful particle definition. Then, (g 2 k (i)) is related to the 
number of produced gravitons and is constant in time. In 
case that amplification of tensor perturbation is still on- 
going, e.g., during a de Sitter phase, the energy density 
related to the massive modes might scale differently The 
scaling behavior (|13[) remains valid also when the fixed 
brane is sent off to infinity and we end up with a single 
braneworld in AdSs, like in the Randall-Sundrum II sce- 
nario 0. The situation is not altered if we replace the 
graviton by a scalar or vector degree of freedom in the 
bulk. Since every bulk degree of freedom must satisfy the 
five-dimensional Klein-Gordon equation, the mode func- 
tions will always be the functions <f) a , and the energy 
density of the KK- modes decays like l/a e . KK particles 
on a brane moving through an AdS bulk cannot play the 
role of dark matter. 

It is important here that we consider a static bulk and 
the time depencence of the brane comes solely from its 
motion through the bulk. In Ref. [l(| the situation of a 
fixed brane in a time- dependent bulk is discussed. There 
it is shown that under certain assumption (separability 
of the y and t dependence of fluctuations), the energy 
density of KK modes on a low energy cosmological brane 
does scale like 1/a 3 which seems to be in contradiction 
with our result. However, the approximations used in [Toj 
lead to a system of equations governing the expansion of 
the Universe but neglecting the time dependence of the 
bulk. The situation is then effectively four dimensional 
even for the KK modes; effects of the fifth dimension like 
the possibility of KK gravitons escaping into the bulk 
seem to be lost in this approach. In our case we would 
have a similar situation if we keep the expansion on the 
brane a(t) but take the position of the brane in the bulk 
as static yb(t) = const, which is not consistent with the 
general relation yb(t) — L/a(t). [For a fixed physical 
mass M = m/a, if we neglect the time dependence of 
<j) n {yb(t)) oc 1/a 2 we also obtain an energy density for 
this mass proportional to 1/a 3 .] 

Particle production: The equation of motion for the 
canonical variables g Q; k is of the form, see Ref. [9j, 

q a ,k + U^klaM = ^ MapqpM + X! Napqp.Vi ■ (14) 

0¥=<* 

Here ui a ,k — \/k 2 + to 2 is the frequency of the mode and 
M. and M are coupling matrices. When we quantize these 
variables, gravitons can be created by two effects: First, 
the time dependence of the effective frequency (w^ ff fe ) 2 = 
w q k~Naa and second, the time dependence of the mode 
couplings described by the antisymmetric matrix M. and 
the off-diagonal part of M . 

Note that Equation (fT4"|) is derived from the corre- 
sponding action for the variables q a ,iL rather than from 
the wave equation ([7]) itself. In this way the approxi- 
mated boundary conditions ((5|) can be implemented con- 



sistently (Km. 

In the technical paper Q we have studied graviton pro- 
duction provoked by a brane moving according to ([3]) in 
great detail numerically. We have found that for long 
wavelengths, kL <C 1, the zero mode is mainly generated 
by its self-coupling, i.e. the time dependence of its ef- 
fective frequency. One actually finds that Aoo <K^)> 
so that there is an instability at the moment of the kink 
which leads to particle creation, and the number of 4D- 
gravitons is given by 2vb/ (kL) 2 . This is specific to ra- 
diation dominated expansion where H 2 a 2 = —d n (Ha). 
For another expansion law we would also obtain particle 
creation during the contraction and expansion phases. 
Light KK gravitons are produced mainly via their cou- 
pling to the zero mode. This behavior changes drastically 
for short wavelengths kL ^> 1. Then the evolution of the 
zero mode couples strongly to the KK modes and pro- 
duction of 4D gravitons via the decay of KK modes takes 
place. In this case the number of produced 4D gravitons 
decays only like oc \/(kL). 

Results and discussion: The numerical simulations 
have revealed a multitude of interesting effects. In the 
following we summarize the main findings. We refer the 
interested reader to Ref. @ for an extensive discussion. 
For the zero-mode power spectrum we find on scales 
kL -C 1 on which we observe cosmological fluctuations 
(Mpc or larger) 

■n i u\ Ki J fc2 if fct « 1 

Po(*) = 2^«6| l {L a)- 2 if kt»l. ^ 

The spectrum of tensor perturbations is blue on super- 
horizon scales as one would expect for an ekpyrotic sce- 
nario. On cosmic microwave background scales the am- 
plitude of perturbations is of the order of (Ho/mp\) 2 and 
hence unobservably small. 

Calculating the energy density of the produced massless 
gravitons one obtains 0] 

n v *> ria\ 
PM ~2?I^- (16) 

Comparing this with the radiation energy density, /9 ra( j = 
(3/(K4-L 2 ))a -4 , the RS fine-tuning condition leads to the 
simple relation 

Pho/prsd — Vb/2. (17) 

The nucleosynthesis bound [12( requests pho ^> 0.1/j ra d, 
which implies Vb < 0.2, justifying our low energy ap- 
proach. The model is not severely constrained by the 
zero-mode. 

More stringent bounds come from the KK modes. Their 
energy density on the brane is found to be 

» 5 «! i 2 

«* - siw (18) 

This result is dominated by high energy KK gravitons 
which are produced due to the kink. It is reasonable 
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to require that the KK-energy density on the brane be 
(much) smaller than the radiation density at all times, 
and in particular, right after the bounce where pkk is 
greatest. If this is not satisfied, back reaction cannot be 
neglected. We obtain with p la d(0) = 3H 2 /k4 



Pkk 

Prad 



a=a(0)=l I \/vb 



Us 



100^ - — . (19) 



L, 



If we use the largest value for the brane velocity Vb ad- 
mitted by the nucleosynthesis bound Vb ~ 0.2 and re- 
quire that PKK/Prad be (much) smaller than one for back- 
reaction effects to be negligible, we obtain the very strin- 
gent condition 



Vs \L 



(20) 



Taking the largest allowed value for L ~ 0.1mm, 
the RS fine-tuning condition Eq. ^ determines 



L s = (LL 2 Pl )^ 3 ~ 10~ 22 mm ~ l/(10 6 TeV) and 
(L/L s ) 2 ~ 10 42 so that y s > L(L/L S ) 2 ~ 10 41 mm 
~ 10 16 Mpc. This is about 12 orders of magnitude larger 
than the present Hubble scale. Also, since yb(t) <C L 
in the low energy regime, and y s 3> L according to the 
inequality (|20|) . the physical brane and the static brane 
need to be far apart at all times otherwise back reaction 
is not negligible. This situation is probably not very 
realistic. We need some high energy, stringy effects to 
provoke the bounce and these may well be relevant only 
when the branes are sufficiently close, i.e. at a distance 
of order L s . But in this case the constraint (|20| will 
be violated which implies that back reaction will be 
relevant. On the other hand, if we want that y s ~ L 
and back reaction to be unimportant, then Eq. (119|) 
implies that the bounce velocity has to be exceedingly 
small, Vb 10~ 15 . One might first hope to find a way 
out of these conclusions by allowing the bounce to 
happen in the high energy regime. But then Vb — 1 and 
the nucleosynthesis bound is violated since too many 



zero-mode gravitons are being produced. Clearly our 
low energy approach looses its justification if Vb — 1, but 
it seems unlikely that modifications coming from the 
high energy regime alleviate the bounds. 
Conclusions: Studying graviton production in an AdS 
braneworld we have found the following. First, the 
energy density of KK gravitons on the brane behaves as 
oc 1/a 6 , i.e. it scales like stiff matter with the expansion 
of the Universe and can therefore not serve as a can- 
didate for dark matter. Furthermore, if gravity looks 
four dimensional on the brane, its higher-dimensional 
aspects, like the KK modes, are repelled from the brane. 
Even if KK gravitons are produced on the brane they 
rapidly escape into the bulk as time evolves, leaving no 
traces of the underlying higher-dimensional nature of 
gravity. This is likely to survive also in other warped 
braneworlds when expansion can be mimicked by brane 
motion. 

Secondly, a braneworld bouncing at low energies is not 
constrained by massless 4D gravitons and satisfies the 
nucleosynthesis bound as long as Vb ^ 0.2. However, 
for interesting values of the string and AdS scales 
and the largest admitted bounce velocity the back 
reaction of the KK modes is only negligible if the 
two branes are far apart from each other at all times, 
which seems rather unrealistic. For a realistic bounce 
the back reaction from KK modes can most likely not 
be neglected. Even if the energy density of the KK 
gravitons on the brane dilutes rapidly after the bounce, 
the corresponding energy density in the bulk could even 
lead to important changes of the bulk geometry. The 
present model seems to be adequate to address the back 
reaction issue since the creation of KK gravitons happens 
exclusively at the bounce. This and the treatment of 
the high energy regime Vb ~ 1 is reserved for future work. 
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